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Abstract 

Forest fires are extremely complex and destructive natural 
phenomena which depend on availability of fuel, 
meteorological conditions and an initial ignition. The overall 
purpose of forest fire research is to make a better 
understanding of the phenomena so that firefighting 
techniques can be improved, controlled burnings can be 
undertaken with more confidence, and lives and property 
might be preserved. In this paper, the assignment and 
theoretical investigations on the problems of crown forest 
fire spread in windy condition were carried out. 
Mathematical model of forest fire was based on an analysis 
of known experimental data and using concept and methods 
from reactive media mechanics. In this context, a study- 
mathematical modeling- of the conditions of forest fire 
spreading that would make it possible to obtain a detailed 
picture of the change in the temperature and component 
concentration fields with time, and determine as well as the 
limiting condition of fire propagation in forest with fire 
break has been carried out. 
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Introduction 

Many mathematical models have been developed to 
calculate forward rates of spread for various fuel 
complexes. A detailed list is given by Grishin 
A.M.(1997). Crown fires are initiated by convective 
and radiative heat transfer from surface fires. 
Fiowever, convection is the main heat transfer 
mechanism (Van Wagner 1977). The theory proposed 
by Van Wagner depends on three simple crown 
properties: crown base height, bulk density of forest 
combustible materials and moisture content of forest 
fuel. Also, crown fire initiation and hazard have been 
studied and modeled in details later by other authors 
(Alexander V.E., 1998;Van Wagner C.E., 1989; 

Xanthopoulos G., 1990; Rothermel R.C., 1991; Cruz 
M.G. et al., 2002; Albini F.A. et al. , 1995; Scott J.H. et 


al., 2001). The more complete discussion of the 
problem of crown forest fires is provided by 
coworkers at Tomsk University (Grishin. 1997, 1998 
and Perminov V.A., 1995). In particular, a 

mathematical model of forest fires was obtained by 
Grishin (1997) based on an analysis of known and 
original experimental data Konev (1977), and using 
concepts and methods from reactive media mechanics. 
The physical two-phase models used by Morvan and 
Dupuy (2001, 2004) may be considered as a 
continuation and extension of the formulation 
proposed by A.M. Grishin(1977). 

This study gives a two dimensional averaged 
mathematical setting and method of numerical 
solution to a problem of a forest fire spread. The 
boundary-value problem is solved numerically using 
the method of splitting according to physical 
processes. It was based on numerical solution of two 
dimensional Reynolds equations for the description of 
turbulent flow taking into account for diffusion 
equations chemical components and equations of 
energy conservation for gaseous and condensed 
phases, volume of fraction of condensed phase (dry 
organic substance, moisture, condensed pyrolysis 
products, mineral part of forest fuel). One aspect of 
this research is to study of the conditions on which the 
forest fire spreads through firebreaks and glades. 

Mathematical Model 

It is assumed that the forest during a forest fire can be 
modeled as 1) a multi-phase, multistoried, spatially 
heterogeneous medium; 2) in the fire zone the forest is 
a porous-dispersed, two-temperature, single-velocity, 
reactive medium; 3) the forest canopy is supposed to 
be non - deformed medium (trunks, large branches, 
small twigs and needles), which affects only the 
magnitude of the force of resistance in the equation of 
conservation of momentum in the gas phase, i.e., the 
medium is assumed to be quasi-solid (almost non- 
deformable during wind gusts); 4) let there be a so- 
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called "ventilated" forest massif, in which the volume 
of fractions of condensed forest fuel phases, consisting 
of dry organic matter, water in liquid state, solid 
pyrolysis products, and ash, can be neglected 
compared to the volume fraction of gas phase 
(components of air and gaseous pyrolysis products); 5) 
the flow has a developed turbulent nature and 
molecular transfer is neglected; 6) gaseous phase 
density doesn't depend on the pressure because of the 
low velocities of the flow in comparison with the 
velocity of the sound. Let the point xi, xi , X3= 0 situated 
at the centre of the surface forest fire source at the 
height of the roughness level, axis Oxi directed parallel 
to the Earth's surface to the right in the direction of the 
unperturbed wind speed, axis 0 x 2 directed 
perpendicular to Oxi and axis 0x3 directed upward (Fig. 
1 ). 



FIG. 1 SCHEME OF THE DOMAIN 
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The system of equations (l)-(6) must be solved taking 
into account the initial and boundary conditions 
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It is considering a constant wind over a perfectly flat 
terrain containing a homogeneous supply of forest 
combustible materials. It was supposed now that there 
is a uniform, constant wind blowing. Because of the 
horizontal sizes of forest massif more than height of 
forest, system of equations of general mathematical 
model of forest fire [1,11] was integrated between the 
limits from height of the roughness level - 0 to /z. 
Besides, suppose that 
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functions for the vertical averaging have been dropped 
in the system of equations, initial and boundary 
conditions for the simplicity of designations. The 
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Here and above — is the symbol of the total 
c it 

(substantial) derivative; a v is the coefficient of phase 

exchange; p - density of gas - dispersed phase, t is 
time; Vi - the velocity components; T, Ts, - 
temperatures of gas and solid phases, Ur - density of 
radiation energy, k - coefficient of radiation 
attenuation, P - pressure; c v - constant pressure specific 
heat of the gas phase, c V i, pi, qx - specific heat, density 
and volume of fraction of condensed phase (1 - dry 
organic substance, 2 - moisture, 3 - condensed 
pyrolysis products, 4 - mineral part of forest fuel), Ri - 
the mass rates of chemical reactions, qx - thermal 
effects of chemical reactions; k g , ks - radiation 
absorption coefficients for gas and condensed phases; 
T e - the ambient temperature; c a - mass concentrations 
of a - component of gas - dispersed medium, index 
cf= 1,2,3, where 1 corresponds to the density of oxygen, 
2 - to carbon monoxide CO, 3 - to carbon dioxide and 
inert components of air; R - universal gas constant; 
M a , Me, and M molecular mass of a -components of 
the gas phase, carbon and air mixture; g is the gravity 
acceleration; Cd is an empirical coefficient of the 
resistance of the vegetation, s is the specific surface of 
the forest fuel in the given forest stratum. In system of 
equations (l)-(6) are introduced by the next 
designations: 

m = pv„T i =-pv' i v' 3 ,J a = -pv[c' a ,J T =-pv[T' 

Upper indexes "+" and designate values of 
functions at X3=h and X3=0 correspondingly. It is 
assumed that heat and mass exchange of fire front and 
boundary layer of atmosphere are governed by 
Newton law and written using the formulas 

(qj -q^)/h = -a(T -T e )/h, 

(J- a -J + a )/h = -a(c-c ae )/hc p . 

To define source terms which characterize inflow 
(outflow of mass) in a volume unit of the gas- 
dispersed phase, the following formulae were used for 
the rate of formulation of the gas-dispersed mixture m , 
outflow of oxygen R 5l , changing carbon monoxide 
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Here v g - mass fraction of gas combustible products of 
pyrolysis, m and as - empirical constants. Reaction 


rates of these various contributions (pyrolysis, 
evaporation, combustion of coke and volatile 
combustible products of pyrolysis) are approximated 
by Arrhenius laws whose parameters (pre-exponential 
constant h and activation energy Ei) are evaluated 
using data for mathematical models (Grishin 1997). 
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The initial values for volume of fractions of condensed 
phases are determined using the expressions: 
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where d -bulk density for surface layer, v z - coefficient 
of ashes of forest fuel, W - forest fuel moisture content. 
It is supposed that the optical properties of a medium 
are independent of radiation wavelength (the 
assumption that the medium is "grey"), and the so- 
called diffusion approximation for radiation flux 
density was used for a mathematical description of 
radiation transport during forest fires. To close the 
system (l)-(6), the components of the tensor of 
turbulent stresses, and the turbulent heat and mass 
fluxes are determined using the local-equilibrium 
model of turbulence (Grishin 1997). The system of 
equations (l)-(6) contains terms associated with 
turbulent diffusion, thermal conduction, and 
convection, and needs to be closed. The components of 
the tensor of turbulent stresses pv-Vj , as well as the 

turbulent fluxes of heat and mass nv'c T , ov’c' are 

r ] p r j a 

written in terms of the gradients of the average flow 
properties using the formulas 
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where Pq, Df are the coefficients of turbulent 

viscosity, thermal conductivity, and diffusion, 
respectively; Pr^, Scf are the turbulent Prandtl and 

Schmidt numbers, which were assumed to be equal to 
1. In dimensional form, the coefficient of dynamic 
turbulent viscosity is determined using local 
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equilibrium model of turbulence (Grishin 1997). The 
thermodynamic, thermophysical and structural 
characteristics correspond to the forest fuels in the 
canopy of a different (for example pine) type of forest. 
The system of equations (l)-(6) must be solved taking 
into account the initial and boundary conditions. The 
thermodynamic, thermophysical and structural 
characteristics correspond to the forest fuels in the 
canopy of a different type of forest; for example, pine 
forest (Grishin and Perminov 1998). 

Numerical Method And Results 

The boundary-value problem (l)-(6) is solved 
numerically using the method of splitting according 
to physical processes. In the first stage, the 
hydrodynamic pattern of flow and distribution of 
scalar functions was calculated. The system of 
ordinary differential equations of chemical kinetics 
obtained as a result of splitting was then integrated. A 
discrete analog was obtained by means of the control 
volume method using the SIMPLE like algorithm (S. 
Patankar, 1981). The accuracy of the program was 
checked by the method of inserted analytical solutions. 
Analytical expressions for the unknown functions 
were substituted in (l)-(6) and the closure of the 
equations was calculated. This was then treated as the 
source in each equation. Next, with the aid of the 
algorithm described above, the values of the functions 
used were inferred with an accuracy of not less than 
1%. The effect of the dimensions of the control 
volumes on the solution was studied by diminishing 
them. The time step was selected automatically. 

Fields of temperature, velocity, component mass 
fractions, and volume fractions of phases were 
obtained numerically. The distribution of basic 
functions shows that the process of crown forest fire 
initiation goes through the next stages. The first stage is 
related to increasing maximum temperature in the fire 
source. At this process stage, the fire source of a 
thermal wind is formed by a zone of heated forest fire 
pyrolysis products which are mixed with air, float up 
and penetrate into the crowns of trees. As a result, 
forest fuels in the tree crowns are heated, as well 
moisture was evaporated, gaseous form and dispersed 
pyrolysis products are generated. Ignition of gaseous 
pyrolysis products of the ground cover occurs at the 
next stage, and that of gaseous pyrolysis products in 
the forest canopy occurs at the last stage. As a result of 
heating of forest fuel elements of crown, moisture 
evaporates, and pyrolysis occurs accompanied by the 
release of gaseous products, which then ignite and 


burn away in the forest canopy. At the moment of 
ignition, the gas combustible products of pyrolysis 
burns away, and the concentration of oxygen is rapidly 
reduced. The temperatures of both phases reach a 
maximum value at the point of ignition. The ignition 
processes is of a gas - phase nature. Note also that the 
transfer of energy from the fire source takes place due 
to radiation; and the value of radiation heat flux 
density is small compared to that of the convective heat 
flux. At Ve ^ 0, the wind field in the forest canopy 
interacts with the gas-jet obstacle that forms from the 
forest fire source and from the ignited forest canopy 
and burns away in the forest canopy. Figures 2-5 
present the distribution of temperature 
f(f -T !T e ,T e = 300/0 (1- 2., 2 - 2.6, 3 - 3, 4 - 3.5, 5 - 
4.) for gas phase, concentrations of oxygen c, (1 - 0.1, 2 
- 0.5, 3 - 0.6, 4 - 0.7, 5 - 0.8, 6 - 0.9) and volatile 
combustible products of pyrolysis c 2 (1-1., 2- 0.1, 3 - 
0.05, 4 - 0.01) (c a = c a / c le , c le = 0.23 ) and temperature 
of condensed phase T s ( T s -T s !T e ,T e = 300AT)(1- 2., 
2 - 2.6, 3 - 3, 4 - 3.5, 5 - 4.) for wind velocity V e = 10 m/s 
at h=10 m: 1) t = 3 sec., 2) t=10 sec, 3) f=18 sec., 4) t= 24 sec. 



FIG. 2 FIELD OF ISOTHERMS OF THE FOREST FIRE SPREAD 
(GAS PHASE). 



FIG. 3 THE DISTRIBUTION OF OXYGEN C x 

It can be noted that the isotherms and lines of equal 
levels are moved in the forest canopy and deformed 
by the action of wind. Similarly, the fields of 
component concentrations are deformed. It is 
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concluded that the forest fire begins to spread. 
Mathematical model and the result of the calculation 
give an opportunity to consider forest fire spread for 
different wind velocity. 



FIG. 4 THE DISTRIBUTION OF C 2 



FIG. 5 FIELD OF ISOTHERMS OF THE FOREST FIRE SPREAD 
(SOLID PHASE) 

Figures 6 ( a , b, c, d) presents the distribution of 
temperature for gas phase, concentration of oxygen 
and volatile combustible products of pyrolysis c 2 
concentrations and temperature of condensed phase 
for wind velocity V e = 5 m/s at h=10 m: 1) t = 3 sec., 2) 

f=10 sec, 3) t = 18 sec., 4) t= 20 sec., 5) t = 31 sec., 6) t= 40 
sec. The results reported in Fig. 6 show the decrease of 
the wind induces a decrease of the rate of fire spread. 
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FIG. 6 FIELDS OF ISOTHERMS OF GAS (A) AND SOLID 
PHASE (D), ISOLINES OF OXYGEN(B) AND GAS PRODUCTS OF 
PYROLYSIS(C) 


One of the objectives of this paper could be to develop 
modeling means to reduce forest fire hazard in forest 
or near towns. In this paper, it presents numerical 
results to study forest fire propagation through 
firebreak. This problem was considered by Zverev 
(1985) in one dimensional mathematical model 
approach. Figures 7 and 8 (Figure 8 b is a continuation 
of Figure 8 a) present the forest fire front movement 
using distributions of temperature at different instants 
of time for two sizes of firebreaks (4.5 and 4 meters). 
The fire break is situated in the middle of domain ( xi 
=100 m). In the first case, the fire could not spread 
through this fire break. 



FIG. 7 FIELD OF ISOTHERMS. FIRE BREAK EQUALS 4.5 M 



FIG. 8 FIELD OF ISOTHERMS. FIRE BREAK EQUALS 4 M 

If the fire break reduces to 4 meters (Figure 8), the fire 
continues to spread but the isotherm (isotherm 5) of 
forest fire is decreased after overcoming fire break. In 
the Figure 9. the dependence of critical fire break value 
for different wind velocities is presented. Of course, the 
size of safe distance depends not only on wind velocity, 
but type and quality of forest combustible materials, its 
moisture, height of trees and others conditions. This 
model allows studying an influence of all these main 
factors. 
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FIG. 9 THE INFLUENCE OF WIND VELOCITY AT THE SIZE 
OF FIRE BREAK. 

Figure 10 (a, b,c) shows the results of numerical 
simulation of a forest fire spreading around the glade 
under the action of wind blowing through it at a speed 
5 m/s in the direction of the Oxi-axis. Initially, the 
source of the fire has the shape of a rectangular. Then 
isotherms are deformed under the action of wind and 
the contour of forest fire looks as crescent (Fig. 10 a, 
curves I). When the fire (isotherms II in Fig. 10 a) moves 
around the forest glade, it is divided into two parts. But 
after that two fire fronts were joined in united fire 
(isotherms III in Fig. 10 a). 



FIG. 10 FIELDS OF ISOTHERMS OF GAS PHASE (A), ISOLINES OF 
OXYGEN(B) AND GAS PRODUCTS OF PYROLYSIS(C) 


Figures 10 (b, c) presents the distribution of 

concentration of oxygen and volatile combustible 
products of pyrolysis c 2 for this case. 

Conclusion 

The results of calculation give an opportunity to 
evaluate critical condition of the forest fire spread, 
which allows applying the given model for preventing 
fires. It overestimates the velocity of crown forest fire 
spread that depends on crown properties: bulk density, 
moisture content of forest fuel and etc. The model 
proposed here gives a detailed picture of the change in 
the temperature and component concentration fields 
with time, and determine as well as the influence of 
different conditions on the crown forest fire spreading 
for the case of inhomogeneous of forest combustible 
materials supply distribution in the area and there 
being such obstacles to fire spread as roads, fire breaks, 
glades, water bodies etc. The results of calculation of 
the rate of crown forest fire spread obtained agree 
with the laws of physics and experimental data for the 
cases when the forest fire spread through the crowns 
of forest under the influence of wind. In this case the 
rate of spread of forest fires is 8-25 km/h (Zverev, 1985 
and Grishin A.M., 1997). 
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